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Abstract.  The  adoption  of  the  International  Celestial  Reference  Sys¬ 
tem  (ICRS)  by  the  IAU  in  use  since  1  January  1998,  and  the  accuracy 
achieved  by  the  most  recent  models  and  observations  of  Earth  rotation 
call  for  a  redefinition  of  the  Earth  Orientation  Parameters  (EOP).  First, 
the  precession-nutation  parameters  and  Greenwich  sidereal  time,  which 
are  currently  defined  in  the  FK5  System,  have  to  be  re-defined  to  be 
consistent  with  the  ICRS.  Second,  the  current  definition  of  the  Celestial 
Ephemeris  Pole  (CEP)  has  to  be  extended  in  order  to  be  consistent  with 
the  most  recent  models  for  nutation  and  polar  motion  at  a  microarsecond 
accuracy  including  diurnal  and  sub-diurnal  components,  as  well  as  with 
new  strategies  of  observations.  Such  issues  have  been  under  considera¬ 
tion  by  the  subgroup  T5  named  “Computational  Consequences”  of  the 
IAU  Working  Group  “ICRS”.  This  paper  gives,  as  the  basis  for  future 
recommendations,  the  preliminary  proposals  of  the  subgroup  T5  for  a 
modern  definition  of  the  CEP,  for  the  definition  of  more  basic  EOP  in 
the  ICRS  and  for  the  choice  of  a  new  origin  on  the  equator  of  the  CEP  in 
place  of  the  equinox.  Then,  the  paper  emphasizes  the  use  of  the  Celestial 
Ephemeris  Origin  (CEO)  which  is  defined  as  the  “non-rotating  origin”  in 
the  celestial  frame  on  the  equator  of  the  CEP. 


1.  Introduction 

There  have  been  important  changes  in  the  situation  of  Earth  rotation  during 
the  last  20  years.  First,  there  have  been  large  improvements  in  the  observations 
since  1980  both  in  precision  and  time  resolution  and  in  strategy  and  procedures. 
An  important  change  is  that  the  celestial  pole  offsets  are  estimated  on  a  regular 
basis  by  VLBI  since  1980  providing  the  actual  position  of  the  pole  in  the  celestial 
frame. 

Second,  the  International  Astronomical  Union  (IAU)  adopted  the  Interna¬ 
tional  Celestial  Reference  System  (ICRS)  in  use  since  1  January  1998  (Ma  et 
al.,  1998)  such  that  the  bary centric  directions  of  distant  extragalactic  objects 
show  no  global  rotation  with  respect  to  these  objects.  The  geocentric  celes¬ 
tial  system  used  for  the  transformation  between  the  celestial  and  the  terrestrial 
system,  designated  here  by  CRS,  is  such  that  its  orientation  follows  the  kine- 
matical  condition  of  absence  of  global  rotation  of  geocentric  directions  of  the 
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objects  that  realize  the  ICRS.  This  provides  a  non-rotating  geocentric  celestial 
frame  to  measure  Earth  rotation  that  has  no  global  rotation  and  is  no  longer 
dependent  on  the  Earth’s  motion  as  was  previously  the  case  with  the  FK5. 

Third,  there  have  been  large  improvements  in  theory  at  pas  accuracy,  dur¬ 
ing  the  period  1990-1999.  The  theory  now  includes  diurnal  and  sub-diurnal 
nutations  in  the  CRS  as  well  as  diurnal  and  sub-diurnal  variations  in  polar 
motion. 

The  consequences  of  this  new  situation  are  such  that  the  current  parameters 
for  precession,  nutation  and  the  Earth’s  angle  of  rotation,  which  are  referred  to 
the  FK5,  are  not  adapted  to  the  ICRS  and  that  the  current  definition  of  the  CEP 
is  not  consistent  with  the  precision  and  spatial  resolution  of  modern  techniques, 
nor  with  the  accuracy  of  theories  and  frequencies  of  the  included  terms.  This 
new  situation  requires  that  more  basic  parameters  for  precession,  nutation  and 
the  Earth’s  angle  of  rotation  be  defined  with  respect  to  the  CRS  and  that  an 
extended  definition  of  the  CEP  be  given. 

Such  consequences  are  under  study  in  the  sub-group  T5  “Computational 
Consequences”  of  the  IAU/WG  ICRS  whose  membership  is:  N.  Capitaine,  V. 
Dehant,  A.-M.  Gontier,  B.  Kolaczek,  J.  Kovalevsky,  J.  Lieske,  C.  Ma,  D.D. 
McCarthy,  0.  Sovers,  and  J.  Vondrak.  There  have  been  four  T5  Newsletters 
from  1998  to  2000,  devoted  to  discussions  of  the  extended  definition  of  the  CEP 
and  of  the  definition  of  new  parameters  for  Earth  orientation.  Other  contributors 
to  the  discussions  of  these  newsletters  have  been:  C.  Bizouard,  P.  Bretagnon, 
A.  Brzezi/’nski,  B.  Guinot,  S.  Loyer,  P.  Mathews,  J.  Mueller,  L.  Petrov,  M. 
Rothacher,  H.  Schuh,  and  Y.  Yatskiv. 

This  paper  gives  the  preliminary  proposals  of  the  subgroup  T5  that  are 
considered  as  the  basis  for  future  recommendations.  The  proposal  for  a  modern 
definition  of  the  CEP  is  discussed  in  Section  2  and  the  definition  of  more  basic 
EOP  in  the  CRS  in  Section  3  together  with  the  choice  of  a  new  origin  on  the 
equator  of  the  CEP.  Then,  the  paper  emphasizes,  in  Section  4,  the  use  of  the 
Celestial  Ephemeris  Origin  (CEO)  which  is  defined  as  the  “non-rotating  origin” 
in  the  CRS  on  the  equator  of  the  CEP,  for  the  definition  of  Earth  rotation  and 
UT1. 

2.  Towards  a  modern  definition  of  the  CEP 
2.1.  Current  realization  of  the  pole 

The  observations  of  the  Earth’s  orientation  are  sensitive  to  the  position  of  the 
pole  of  the  Terrestrial  Reference  System  (TRS)  with  respect  to  the  CRS.  How¬ 
ever,  an  intermediate  pole  is  considered  in  the  coordinate  transformation  be¬ 
tween  the  TRS  and  the  CRS.  This  is  currently  the  Celestial  Ephemeris  Pole 
(CEP)  as  defined  in  the  IAU  1980  theory  of  nutation  (Seidelmann,  1982).  The 
current  Earth’s  angle  of  rotation  as  used  for  estimating  UT1  is  about  the  axis 
of  the  CEP. 

The  CEP  is  an  intermediary  pole  separating,  by  convention,  the  motion  of 
the  pole  of  the  TRS  (i.e.  the  pole  of  the  mean  surface  geographic  figure  axis,  or 
Tisserand  mean  axis  of  the  Earth)  into  two  parts: 
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-  a  “celestial  part”,  designated  as  the  celestial  motion  of  the  CEP,  including 
all  the  long-period  terms  in  the  CRS  (precession/nutation)  with  periods 
greater  than  2  days  (he.  frequencies  between  —1/2  cpd  and  +1/2  cpd), 

-  a  “terrestrial  part”,  designated  as  the  terrestrial  motion  of  the  CEP,  in¬ 
cluding  all  the  long-period  terms  in  the  TRS  (polar  motion)  with  periods 
greater  than  2  days,  (he.  frequencies  between  —1/2  cpd  and  +1/2  cpd). 

Such  motions  can  be  viewed  either  in  the  CRS  or  in  the  TRS  according  to 
the  chosen  convention  as  shown  in  the  schematic  representation  below. 

as  viewed  in  the  TRS _ | _ | _ | _ | _ | _ 

-3/2  -1  -1/2  0  +1/2  frequency  (cpd) 

Celestial  motion  I  Terrestrial  motion 


as  viewed  in  the  CRS _ | _ | _ | _ | _ | _ 

-1/2  0  +1/2  1  +3/2  frequency  (cpd) 


It  must  be  noted  that  there  is  no  convention  for  terms  outside  the  fre¬ 
quency  interval  [—3/2,  +1/2]  as  viewed  in  the  TRS  or  [—1/2,  +3/2]  as  viewed 
in  the  CRS.  It  means  that  the  current  definition  of  the  CEP  is  valid  only  in  this 
frequency  interval  and  does  not  include  any  convention  for  the  high-frequency 
motions  in  both  the  CRS  and  the  TRS.  Moreover,  the  current  convention  is  only 
valid  for  sampling  intervals  greater  than  1  day. 

2.2.  The  need  for  an  extended  definition  of  the  CEP 

The  main  changes  for  the  definition  of  the  CEP  (see  Capitaine  &  Brzezinski, 
1999  or  Capitaine,  1999  for  more  details)  concern  the  predictable  high-frequency 
components  in  nutation  (Bretagnon  et  al. ,  1997)  and  in  polar  motion  (Herring  & 
Dong,  1994)  and  the  new  strategy  of  observations  in  which  sub-daily  observations 
will  increase  in  the  near  future.  The  recent  models,  at  a  microarsecond  accuracy, 
include: 

-  prograde  semidiurnal  nutations  with  amplitudes  of  the  order  of  15  mi- 
croarseconds,  produced  by  the  lunisolar  torque  exerted  on  the  tesseral 
coefficients  (2,2)  and  (3,2)  of  the  Earth’s  gravitational  potential, 

-  prograde  diurnal  nutations  (in  Ar/isineo  and  Ae)  with  amplitudes  of  the 
order  of  15  microar seconds,  due  to  the  tesseral  coefficients  (3,1)  and  (4,1), 

-  prograde  diurnal  tidal  variations  in  polar  motion  with  amplitudes  of  the 
order  of  150  microar  seconds, 

-  prograde  semidiurnal  tidal  variations  in  polar  motion  with  amplitudes  of 
the  order  of  100  microar  seconds, 

-  retrograde  diurnal  tidal  variations  in  polar  motion  with  amplitudes  of  the 
order  of  800  microar  seconds, 
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-  retrograde  semidiurnal  tidal  variations  in  polar  motion  with  amplitudes  of 
the  order  of  300  microar seconds. 

The  overlapping  between  these  high  frequency  motions  is  such  that: 

-  the  prograde  diurnal  nutations  in  space  appear  in  the  Earth  as  long  peri¬ 
odic  prograde  and  retrograde  variations  in  polar  motion  of  periods  ±  27 
d,  193.6  d,  3193.9  d,  ...  (Bizouard  et  al,  1999a), 

-  the  semidiurnal  prograde  nutations  in  space  appear  as  prograde  diurnal 
variations  in  polar  motion,  also  called  “polar  libration”  (Chao  et  al .,  1991), 

-  the  retrograde  diurnal  terms  in  the  tidal  polar  motion  variations  are  in 
fact  included  in  the  most  recent  models  of  nutation  for  a  nonrigid  Earth. 

The  extended  definition  of  the  CEP  must  therefore  clearly  specify  which 
part  of  these  high  frequency  components  is  considered  in  the  CRS  or  in  the  TRS 
on  the  basis  of  clear  theoretical  considerations.  For  sub-daily  observations,  the 
two  frequency  intervals  of  polar  motion  and  nutation  are  no  longer  separated 
and  thus,  the  frequency  domain  decomposition  into  polar  motion  and  nutation 
becomes  strongly  dependent  on  the  method  used  for  processing  the  observations 
(Brzezihski  &  Capitaine,  1993).  Hence,  simultaneous  estimation  of  polar  motion 
and  nutation  must  be  avoided. 

2.3.  Proposal  for  a  modern  definition  of  the  CEP 

The  discussions  in  T5  have  led  to  the  following  proposals  concerning  the  concept 
and  realization  of  the  CEP: 

Concept 

-  to  abandon  the  IAU  1980  definition  “no  nearly- diurnal  motion...”, 

-  the  pole  must  not  be  defined  by  its  realization  but  a  clear  concept  not 
dependent  on  further  improvements  in  the  model, 

-  the  pole  must  be  defined  such  that  it  can  be  realized  by  a  model  as  accu¬ 
rately  as  necessary, 

-  the  change  from  the  current  CEP  should  have  minimal  impact  on  users, 

-  to  consider  a  change  of  name  and  call  the  new  pole  “Celestial  Reference 
Pole”  (CRP),  or  “Celestial  Intermediate  Pole”  (CIP), 

-  to  specify  the  motion  in  the  CRS  by  a  model, 

-  to  sharpen  the  definition  of  the  CEP  by  taking  into  account  the  high 
frequency  variations  of  polar  motion. 

A  tentative  definition  is  the  Pole  of  the  intermediate  equator  (between  the  TRS 
and  CRS )  the  motion  of  which,  in  the  CRS,  is  produced  by  the  external  gravita¬ 
tional  forces  acting  on  the  Earth  after  filtering  out  terms  with  periods  less  than 
2  days. 
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Figure  1.  Precession-nutation  and  Earth’s  rotation  parameters  re¬ 
ferred  to  the  FK5. 


The  forced  motion  of  the  CEP  in  the  C'RS  is  that  of  the  Tisserand  mean 
axis  of  the  Earth.  The  nutations  with  periods  less  than  2  days  (prograde  diurnal 
and  prograde  semidiurnal  nutations)  are  included  in  the  model  for  the  variations 
in  polar  motion. 

Realization 

-  the  motion  of  the  CEP  in  the  C'RS  is  specified  by  the  best  available 
precession-nutation  series  plus  the  time  dependent  celestial  pole  offsets, 

-  the  motion  of  the  CEP  in  the  TRS  is  provided  by  observations  taking  into 
account  a  predictable  part  specified  by  an  empirical  model, 

-  the  corrections  to  the  empirical  model  can  be  obtained  by  extracting  the 
high  frequency  signal  from  the  pole  coordinates  either  together  with  the 
long  periodic  motion  (following  the  procedure  proposed  by  Mathews  ( 1999) 
and  successfully  tested  by  Bizouard  et  al.  (1999b)),  or  in  a  second  step. 


3.  The  Earth  Orientation  Parameters 
3.1.  The  current  EOP 

The  FK5,  based  on  positions  and  proper  motions  of  bright  stars,  was  oriented 
so  that  at  the  “epoch”,  the  positions  are  referred  to  the  best  estimate  of  the 
location  of  the  mean  pole  and  mean  equinox.  The  proper  motions  of  stars  were 
evaluated  so  that,  for  the  adopted  model  of  precession,  they  provide  the  best 
access  to  the  mean  pole  and  mean  equinox  at  any  other  date. 

The  current  precession  angles  (see  Fig.  1)  are  those  defined  by  Lieske  et  al. 
(1977)  in  the  FK5  system.  The  current  nutation  angles  At/’  in  longitude  and  Ac 
in  obliquity  are  referred  to  the  ecliptic  of  date  (see  Fig.  1)  and  \n  +  A\  is  the 
angular  distance  between  the  ecliptic  of  epoch  and  the  ecliptic  of  date  along  the 
equator  of  date. 

The  current  procedure  to  account  for  precession  and  nutation  is  to  use  the 
matrix  transformation  P(t)  using  the  developments  as  function  of  time  of  the 
precession  angles,  C,a-,Qa-,za,  followed  by  the  matrix  transformation  N(t)  using 
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the  nutation  angles  A  ip  and  Ac  provided  by  the  conventional  series  of  nutation 
(Seidelmann,  1982;  McCarthy,  1996).  Such  a  transformation  corresponds  to 
a  sequence  of  six  consecutive  rotations  for  precession  and  nutation  using  five 
different  parameters.  As  the  precession  and  nutation  angles  are  referred  to 
the  ecliptic  of  date,  the  PN(t)  matrix  is  dependent  both  on  the  precession  and 
nutation  of  the  equator  and  to  the  precession  of  the  ecliptic. 

The  current  procedure  to  account  for  Earth  rotation  in  the  FK5  system  is  to 
use  the  relationship  between  Greenwich  sidereal  time  and  UT1  (Aoki  et  al. ,  1982) 
giving  GMST  at  date  /,  followed  by  the  relationship  between  GST  and  GMST 
and  then  to  take  into  account  the  interval  of  GMST  from  Oh  UT1  to  the  hour  of 
observation.  The  difference  GST  — GMST  is  provided,  since  1  January  1997,  by 
the  “complete”  equation  of  the  equinoxes  (Aoki  &  Kinoshita,  1983;  McCarthy, 
1996).  In  addition  to  Earth  rotation,  the  angle  GST  thus  includes  (Capitaine 
&  Gontier,  1993)  a  part  due  to  the  accumulated  precession  and  nutation  along 
the  equator  as  well  as  a  part  (in  GST  — GMST)  due  to  crossed  terms  between 
precession  and  nutation  and  crossed  nutation  terms.  It  refers  to  the  ecliptic  of 
date  and  thus  mixes  Earth  rotation  and  precession-nutation.  It  must  also  be 
noted  that  the  secular  term  in  the  relationship  between  GMST  and  UT1  mixes 
UT1  and  TCG. 

3.2.  Earth  Orientation  Parameters  referred  to  the  CRS 

The  fundamental  properties  of  the  CRS  are  the  absence  of  global  rotation  and 
the  abandonment  of  the  link  with  the  motion  of  the  Earth.  In  this  spirit,  the 
definition  of  the  EOP  must  be  based  on  location  and  motion  of  the  rotation 
pole  only,  without  any  relation  to  the  orbital  motion  of  the  Earth.  It  is  thus 
necessary  to  abandon  the  current  parameters  in  the  FK5  System  which  com¬ 
bine  the  motions  of  the  equator  and  the  ecliptic  with  respect  to  the  ICRS.  The 
requirements  for  the  new  EOP  are  the  following: 

-  to  be  based  on  clear  concepts  referred  to  the  ICRS, 

-  to  use  a  minimal  number  of  parameters, 

-  not  to  be  referred  to  the  ecliptic  of  date,  but  to  a  fixed  plane, 

-  to  abandon  the  use  of  the  true  equinox  as  the  origin  on  the  moving  equator, 

-  to  include  both  precession  and  nutation  of  the  equator. 

Several  options  can  been  considered  among  which  the  most  basic  ones,  the 
Euler  angles  between  the  CRS  and  the  TRS  and  the  celestial  and  terrestrial 
coordinates  of  the  CEP,  are  described  below. 

(i)  Euler  angles  between  the  CRS  and  the  TRS 

The  two  first  Euler  angles,  0  and  ip,  include  both  polar  motion  and  precession- 
nutation  referred  to  the  ecliptic  of  epoch  and  the  third  one,  <p>,  the  origin  of 
which  is  the  node  of  the  equator  of  figure  on  the  ecliptic  of  epoch,  includes  both 
Earth  rotation  and  the  precession-nutation  motion  of  the  equator  of  figure.  No 
intermediate  pole  is  used  in  the  transformation  between  the  TRS  and  the  CRS 
using  Euler  angles  which  is  based  on  only  three  parameters. 
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Figure  2.  The  Euler  angles  between  the  terrestrial  frame  (OABC) 
and  the  celestial  ecliptic  frame  (OXYZ)  (Woolard,  1953). 


z„ 


Figure  3.  The  coordinates  of  the  CEP  (noted  P)  in  the  C'RS  and  TRS. 
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Figure  4.  Alternative  origins  on  the  moving  equator. 


(ii)  The  coordinates  of  the  CEP 

The  coordinates  A”  =  sin  d  cos  P,  Y  =  sine?  sin  A  of  the  CEP  in  the  C'RS  (Fig. 
3)  include  precession  and  nutation,  the  coupling  effects  between  precession  and 
nutation,  as  well  as  the  offsets  £0,  ?/0,  between  the  CEP  at  J2000.0  and  the  pole 
of  the  C'RS  and  dap  between  the  mean  equinox  at  J2000.0  and  the  origin  on  the 
equator  of  the  C'RS.  A”  and  Y  can  be  expressed  as  (C'apitaine,  1990): 

A”  =  A  +  £0  —  dcx0Y ,  Y  =  Y  +  ?/0  +  dap  A”,  with  (see  Fig.  1) 

A  =  sin  ( cjp  +  Auq )  sin(  r/’q  +  Ar/q ) 

Y  =  —  sin  e0  cos(ax4  +  Auq  )  +  cos  e0  siiffaxr  +  Auq  )  cos(  r/’q  +  At|q ), 
the  quantities  r/’q  +  Ar/q  and  0x4  +  Auq  being  referred  to  the  ecliptic  of  epoch. 

The  celestial  pole  coordinates  are  supposed  to  be  referred  to  an  intermediate 
pole,  such  as  the  CEP,  and  appear  in  a  symmetric  form  as  the  coordinates  in 
the  TR.S,  u  =  sin  g  cos  P,  v  =  sin  g  sin  F. 

(iii)  The  choice  of  an  origin  on  the  moving  equator  The  possible  origins  on 
the  moving  equator  to  be  used  in  place  of  the  equinox  7  (intersection  with  the 
moving  ecliptic Jijqare  (see  Fig.  4  and  Fig.  5): 

:  the  intersection  with  the  fixed  ecliptic, 

Iv  :  the  intersection  with  the  zero-line  meridian  of  the  C'RS, 

H  :  the  intersection  with  the  instantaneous  prime  meridian  PSo, 

S  :  such  that  :  SfV  =  So IV, 

N  :  the  node  upon  the  celestial  reference  great  circle, 

a  :  the  non  rotating  origin  (Guinot,  1979),  such  that  i  =  E(cosd  —  1). 

A  comparison  of  the  parameters  considered  above  shows  that  the  Euler 
angles,  or  equivalently  the  celestial  coordinates  of  the  z-axis  of  the  TR.S,  which 
reduce  to  three  the  number  of  EOP,  do  not  use  any  intermediate  pole.  The 
Euler  angles  thus  include  both  high-frequency  and  low-frequency  components 
of  the  motion  of  the  z-axis  of  the  TR.S  in  the  C'RS,  arid  their  estimation  from 
observations  would  need  a.  large  number  of  parameters. 


Definition  of  the  CEP  and  the  CEO 


161 


Figure  5.  Definition  of  the  CEO,  <7,  the  TEO,  w  and  the  stellar 
angle  0. 


On  the  contrary,  the  celestial  coordinates  of  the  CEP  separate  the  celestial 
components  from  the  terrestrial  components  according  to  a  “frequency  criterion” 
using  an  intermediate  pole.  Such  a  procedure  facilitates  the  estimation  of  the 
parameters  from  observations. 

A  comparison  of  the  origins  considered  above,  either  of  geometrical  meaning 
(such  as  S  or  K)  or  of  kinematical  origin  (o)  shows  that,  in  the  first  case,  the 
instantaneous  rotation  of  the  origin  considered  is  included  in  the  derived  value 
for  the  Earth’s  angular  rotation,  whereas  it  is  clearly  separated  when  using 
directly  the  non-rotating  origin  (NRO)  as  the  origin  on  the  moving  equator. 
The  preferred  parameters  resulting  from  the  discussions  of  T5  are  the  celestial 
pole  coordinates  X,  Y,  and  the  preferred  origin  is  the  NRO.  Section  4  emphasizes 
the  practical  use  of  the  NRO  in  the  transformation  coordinates  between  the  TRS 
and  the  C'RS. 


4.  Definition  and  use  of  the  celestial  and  terrestrial  ephemeris  origins 
4.1.  Definition  of  the  CEO  and  TEO 

To  define  an  Earth’s  angle  of  rotation  that  includes  only  the  “intrinsic  Earth 
rotation”,  it  is  necessary  to  eliminate  the  spurious  rotation  of  the  origin  on 
the  moving  equator  due  to  the  precession-nutation  of  the  equator.  It  is  thus 
necessary  to  use  an  origin  defined  by  the  kinematical  property  of  having  “no 
instantaneous  rotation”  around  the  axis  of  rotation  with  respect  to  the  C'RS. 
That  provides  the  “non-rotating  origin”  a  (Guinot,  1979),  which  can  be  des¬ 
ignated  as  “the  Celestial  Ephemeris  Origin”  (CEO)  when  referred  to  the  axis 
of  the  CEP.  Similarly  the  “Terrestrial  Ephemeris  Origin”  (TEO),  m,  designates 
the  terrestrial  NRO.  The  “stellar  angle”  is  defined  by  9  =  zua  (see  Fig.  5)  and 
is  such  that  it  includes  only  Earth  rotation. 

The  positioning  of  the  CEO  and  TEO  are  provided  by  the  quantities  s  and 
s'  such  that,  with  <r(/o)No  =  SoNo: 

s  =  <tN  —  S0N  =  f  ( cos  d  —  1)E  dt  =  —  f  X\+\X  dt  , 

to  to 

t  t 

s'  =  zaM  -  n0M  =  -  faz^-dt  =  f  ( cos g  —  1)F dt  . 

to  to 
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4.2.  The  use  of  the  CEO  and  the  TEO 

The  coordinate  transformation  from  the  TRS  to  the  CRS  using  the  celestial  pole 
coordinates,  the  CEO  and  TEO,  is  composed  of  three  components  W{t),  R(t ) 
and  PN(t ),  respectively  for  polar  motion,  rotation  of  the  Earth  and  precession- 
nutation.  They  can  be  written  as: 

W(t)  =  R3(-sr)  ■  R1(-v)  ■  R2(u)  , 

R(t)  =  R3{-9)  , 

PN(t)  =  R3(-E)  •  R2(-d)  •  R3(E)  •  R3(,s)  . 

Development  as  a  function  of  time  of  X  and  Y  and  s  with  microarsecond  ac¬ 
curacy  after  one  century  are  provided  by  Capitaine  et  al.  (2000)  consistent  with 
the  IERS  1996  precession  and  nutation  and  the  numerical  values  of  the  celestial 
offsets  at  J2000.0.  Such  developments  include  a  constant  term,  a  polynomial 
form  of  /,  a  sum  of  periodic  terms  and  a  sum  of  Poisson  terms.  Using  such 
parameters,  the  definition  of  UT1  takes  the  simple  form:  UT 1  —  UT lo  =  k6 , 
where  k  is  a  constant  selected  so  that,  on  the  average,  the  Sun  crosses  the  prime 
meridian  at  12  hours  UT1.  This  shows  that  the  use  of  the  CEO  and  the  TEO 
provides  EOP  corresponding  to  clear  concepts  referred  to  the  ICRS. 
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